Astronomy & Astrophysics manuscript no. 


February 2, 2008 


(DOI: will be inserted by hand later) 





The kinematics of the bi-lobal supernova remnant G 65.3+5.7 - 

Paper II 

p. Boumis^, J. Meaburn^, J. A. Lopez^, F. Mavromatakis^ , M. P. Redman^, D.J. Harman^, and C. D. 

Goudisi'6 

' ^ Institute of Astronomy & Astrophysics, National Observatory of Athens, I. Metaxa & V. Paulou, P. Penteli, 

^ ■ GR- 15236 Athens, Greece. 

^ Instituto de Astronomia, UNAM, Apdo. Postal 877. Ensenada, B.C. 22800, Mexico. 
^ ' * University of Crete, Physics Department, P.O. Box 2208, 710 03 Heraklion, Crete, Greece. 

^ Dublin Institute for Advanced Studies, School of Cosmic Physics, 5 Merrion Square, Dublin 2, Republic of 
Ireland. 

^ Astronomical Laboratory, Department of Physics, University of Patras, 26500 Rio-Patras, Greece. 
Received 9 March 2004 / Accepted 11 May 2004 

> 

' Abstract. Further deep, narrow-band images in the light of [O ill] 5007 A have been added to the previous mosaic 

l/^ , of the faint galactic supernova remnant G 65.3+5.7. Additionally longslit spatially resolved [O ill] 5007 A line 

CO ' profiles have been obtained at sample positions using the Manchester Echelle Spectrometer at the San Pedro 

I Martir observatory. The remnant is shown to be predominantly bi-lobal with an EW axis for this structure. 

However, a faint additional northern lobe has now been revealed. 

Splitting of the profiles along the slit lengths, when extrapolated to the remnant's centre, although uncertain 
suggests that the expansion velocity of this remnant is between 124 and 187 km s~^ ie much lower than the 400 
■ km s~^ previously predicted for the forward shock velocity from the X-ray emission. 

An expansion proper motion measurement of 2.1 it 0.4 arcsec in 48 years for the remnant's filamentary edge in 
the hght of Ha + [N ll] 6548 & 6584 A has also been made. When combined with an expansion velocity of « 155 
5_j ' km s~^, a distance of « 800 pc to G 65.3+5.7 is derived. 

^ ' Several possibilities are considered for the large difference in the expansion velocity measured here and the 400 

. km s~^ shock velocity required to generate the X-ray emission. It is also suggested that the morphology of the 

remnant may be created by a tilt in the galactic magnetic field in this vicinity. 
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1. Introduction clumpy X-ray emission from G 65.3+5.7, assuming a dis- 
tance of 1 kpc, from a forward shock of « 400 km s^^ in 

Mavromatakis et al. (2002 - hereafter called Paper I) pre- ^ ^^^^^^^ ^^^-^^^ ^^^^^ ^^^^.^^ 02 cm-^. 

sented the deepest optical images of the supernova rem- 

nant (SNR) G 65.3+5.7, discovered by Gull et al. (1977), Losinskaya (1980), from Ha interferograms obtained 

yet obtained. They revealed a filamentary structure « 3° ^^^^ ^ ^^^^^^^^^ Fabry-Perot interferometer, had extrap- 

X 4° in size, emitting the [O III] 5007 A nebular hue ^''gfg^'gVSir'rnewf^^^^^^^ 

brightly and predominantly composed of two irregular fil- ^ th ' t~'~^''^ (called newly discovered SNR in Cygnus 
amentary rings with displaced centres; a bi-lobal morphol- ^ paper j o give a g o a expansion ve oci y cxp — 

. . 5 400 + 200 km s^^. 
ogy IS implied. 

Also presented in Paper I are low-dispersion spectra Deep, spatially resolved profiles of the [O ill] 5007 A 
at eight positions (Pl-8) around the remnant's perimeter; emission line have now been obtained with an echelle spec- 
these are consistent with ionization by a 90-140 km s'^ tromctcr at positions Pl-5 & 7-8 where the Paper I low- 
shock with a local ionized gas density of 200 cm'^ typi- dispersion spectra were obtained, and at a further posi- 
cal of an « 20,000 yr old SNR. ROSAT observations by ^0^' P9' o^er the centre of G 65.3+5.7 with the inten- 
Aschenbach (1994) and Lu & Aschenbach (2004) reveal ti^^ of measuring the maximum expansion velocity di- 
rectly as well as obtaining accurate line profiles at repre- 
hend offprint requests to: Dr. P. Boumis sentive positions. Further deep CCD images, again in the 
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[O III] 5007 A line and neighbouring continuum, were also 
obtained, with the same wide-field imaging system em- 
ployed in Paper I; a more complete continuum subtracted 
image of this SNR has resulted. 

2. Observations and results 

2.1. New imagery 

Three further [O III] 5007 A and continuum images 
were obtained on June 2003 27-28 with the 89' x 89' 
field (5" per pixel) 0.3 m Schmidt Cassegrain telescope 
at the Skinakas Observatory, Crete, Greece to complete 
the coverage of the westerly edge of the SNR. These 
were then combined, after an identical analysis, with the 
[O III] 5007 A image described in Paper I. Light and deep 
negative grey-scale representations of the whole field are 
shown in Figs, la-b respectively. 

2.2. Long-slit spectroscopy 

The present spectral observations were made with 
the Manchester Echelle Spectrometer (MES-SPM - see 
Meaburn et al. 1984; 2003) combined with the 2.1-m San 
Pedro Martir telescope. A SITe3 CCD was the detector 
with 1024x1024, 24 /xm pixels although 2x2 binning was 
employed throughout the observations on the nights of the 
May 2003 29-31. 

Spatially resolved, long-slit line profiles at high spec- 
tral resolution were obtained with the MES-SPM. This 
spectrometer has no cross-dispersion. For the present ob- 
servations, a filter of 60 A bandwidth was used to isolate 
the 114"^ echelle order containing the [O iii] 5007 A neb- 
ular emission line. 

The 512 increments, each 0.626" long, give a total 
projected slit length of 5^34 on the sky. 'Seeing' was al- 
ways < 1" during these observations. A 150 fim wide 
(= 12 kms~^ and 1.9") single slit was used. 

The data were bias corrected, cleaned etc. in the usual 
way using the STARLINK FIGARO and kappa software 
packages. All spectra were calibrated in heliocentric radial 
velocity ( Fhel) to ± 3 kms~^ accuracy against spectra 
of a thorium/argon lamp. Absolute surface brightnesses, 
B[o III] erg s~^ cm~^ sr~^ A~^, of the line profiles were 
obtained by comparing the spectra to the slitless spectrum 
of the standard star Feige 56. All spectra were obtained 
in photometric conditions, and without correction for in- 
terstellar extinction are accurate to around ±10 percent. 

The slit was orientated EW on the centres Pl-5 & 
P7-8 from Paper I with integration times of 1800s. The 
integration time for a further slit position (P9) was 3600s. 
These slit positions are shown against the lightly printed 
image of the SNR in Fig. la. 

Negative grey-scale representations of the position- 
velocity (pv) arrays of [O ill] 5007 A line profiles are shown 
in Figs 2a h for slit positions Pl-5 & P7 8. More details 
of the imagery and pv arrays over the brighter filaments 



are shown in Figs. 3a-b and 4a-b for slits PI and P2 re- 
spectively. 

The line profiles extracted from various lengths of the 
pv arrays in Fig. 2a-h and for P9, where the emission is 
very faint, are shown in Fig. 5. The observed peak value 
of B[o III] in units of 10~^ erg s~"^ cm~^ sr"-"^ A~-^ is given 
next to each profile. 

2.3. Expansion Proper Motions 

The filament at Pos. 2 (Fig. la) was detected by the 
Paloniar Observatory Sky Survey POSS-E red plate taken 
in 1951. This is 48 yr before the Ha + [N ii] 6548 & 6584 A 
image was obtained (taken in 1999 July 13 similarly to 
those in Fig. la & b - see Paper I for details). Both 
images are predominantly in the light of the Ha -1- 
[N ii] 6548 & 6584 A lines and the Digitized Sky Survey 
of the POSS plate has a very similar angular resolution to 
the 1999 imagery. A measurement of the expansion proper 
motion of G 65.3-1-5.7 has therefore been made. 

Firstly, the two time-separated images were rotated 
identically until the filament covered by Pos. 2 became 
vertical in both arrays. Profiles were then extracted along 
identical lines where there were no confusing stellar im- 
ages on the filament but many reference star profiles fur- 
ther away. Gaussians were fitted to both faint stellar and 
filament profiles to give a measured shift perpendicular 
to the filament's length, and away from the centre of the 
supernova remnant, of 69pm = 2.1 ± 0.4 arcsec in 48 years. 

3. Discussion 

3.1. Kinematics and distance 

The first useful parameter to determine for such an ex- 
tended SNR is its systemic heliocentric radial velocity 
( ^ys)- This is given most plausibly by the mean Vhel 
of the profiles of the brightest filaments on the perime- 
ter of the remnant for it is assumed that these are being 
viewed tangentially. The centroids of the single profiles in 
Figs. 3, 4 & 5 for the filamentary edges along slits PI, 
2, 3, 5 and 7 are at Vhel = -19, 0, 0, -10 and -15 km 
s^^ respectively to give a value of Vsys = -7 km s^^ for 
G 65.3+5.7 (i.e. halfway between the mean radial veloci- 
ties of the northern and southern filaments). Note that the 
profiles for the northern filaments (PI, P7 and P5) are sig- 
nificantly displaced to approaching radial velocities com- 
pared with their southern counterparts (P2 & P3). The 
halfwidths of the same profiles (for PI, 2, 3, 5 & 7) fit- 
ted by single Gaussians and corrected for the instrumental 
broadening are 0.42 ± 0.02, 0.44 ± 0.03, 0.47 ± 0.01, 0.46 
± 0.08 & 0.36 ± 0.01 A respectively As the post-shock 
interstellar gas emitting [O III] 5007 A has most likely 
cooled to Te = 10^ K then these widths imply residual 
turbulent motions, c;ombined with non-turbulent motions 
due to curvature of the shock fronts, of around 23 km 
s~^ in all of these filamentary edges for the width of the 
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thermal component of the [O ill] 5007 A hne wiU only be 
0.09 A as given by 6X = 8.927 x 10"^ x T^^ A. 

Knowledge of the global expansion velocity Vcxp for 
such a complex remnant is also desirable. In the pv ar- 
rays in Fig. 2, as detected by Lozinskaya (1980), there is 
evidence of expansion along slits P2 and P7: the profiles 
are split by about 70 km s~^ away from the bright fil- 
amentary edge, towards the westerly end of P2 (and see 
this trend along the short length shown in Fig. 4b) and 
go from -15 to -60 km s^^ along the length of P7. If it 
is assumed that this splitting towards the westerly end 
of P2 is due to the spherical expansion of the 91' radius 
easterly lobe of G 65.3-1-5.7 then, after taking account of 
the slit orientation with respect to the filamentary edge, 
by extrapolation Vcxp = 187 ± 50 km is derived for 
this lobe of the remnant. The lower end of the Kxp values 
given by Losinskaya (1980) ( to ± 200 km s~^ accuracy) 
is therefore favoured. 

Other clues to the true value of V^xp arc also in the 
present observations. For instance velocity components 
further from T^ys are found along slit positions P8 and 9 
which are well inside the remnant's outer perimeter (Fig. 
la). In the pv array in Fig. 2 for P8 a component at Vhel 
» 75 km s^^ can be seen. If this reflects the expansion of 
the easterly lobe in Fig. la, and assuming Vsys = -7 km 
s""*^ and spherical expansion, then V^^p ~ 124 km s~^ 
for the lobe for this position is w 0.25 of the remnant's 
radius towards the lobe's centre; whereas the component 
with Vhel = km s^^ along the same slit length must be 
from the filamentary edge of the western lobe (see Fig. 1 
a). Similarly, the faint profile from P9 in Fig. 5 from near 
the center of the eastern lobe but just inside the edge of 
the western lobe of the remnant has a velocity component 
at Vhel = +130 km s^^ which when compared with Vgys 
= -7 km s^^ could imply V^xp ~ 140 km s^^ maybe for 
the eastern lobe and most likely for the whole remnant. 

All of the current estimations of V^xp (between 124 to 
187 km s~^) are well short of the value of the 400 km s~^ 
forward shock velocity required for the X-ray emission (Lu 
& Aschenbach 2004) and below the bottom end ( V^xp = 
200 km s~^) of the large range given by Losinskaya (1980). 
If spherical expansion is assumed at V^xpi km s~^) in T 
(yr) to give 56pm (arcscc) then the distance D(pc) is given 
by D = 0.2168 x Vexp x T x 69 in which case for d9pm 
= 2.1 arcsec in 48 yr then D/ Vexp = 4.955. This gives D 
= 770 ± 200 pc for Fexp » 155 km s'^ but a large 1980 
pc for 400 km s~^. 

3.2. Energetics 

As discussed above the optically measured expansion and 
shock velocities and the X-ray derived shock velocities dif- 
fer sharply. This is also observed for other similar SNR, 
particularly the Cygnus Loop. Lu & Aschenbach (2004) 
and others have assumed that G65.3-I-5.7 is still in its adi- 
abatic phase of evolution, thus enabling a Sedov analysis of 
the energetics. The low ambient density derived from this 



analysis suggests a lower than usual post-shock cooling 
rate (oc n?) and thus make it less than likely that the SNR 
will have moved into the momentum conserving phase. 
The simplest solution to the discrepancy is that the cen- 
tral regions were indeed subject to a ~ 400 km s~^ shock 
in a rarified medium which produced the X-ray emission 
but that the shock speed has recently encountered much 
denser material which has reduced the shock speed and 
led to the optical emission that is observed. Two possi- 
ble sources of this dense material are a cloudy ISM and 
a pre-existing dense circumstellar shell and these are now 
discussed in turn. 

Lu & Aschenbach (2004) favour a shock wave propa- 
gating into a cloudy ISM as described in detail by McKee 
& Cowie (1975). The two shock speeds are then due to 
cloud shocks and the main shock travelling through the 
intercloud medium. Less obviously explained is the fact 
that the measured expansion velocity of the SNR as mea- 
sured by the optical filaments is also much less than the 
X-ray shock velocity. The X-ray emission does not extend 
measurably beyond the boundary delimited by the opti- 
cal filaments. This requires that the clouds are promptly 
disrupted by the shock, incorporated into the postshock 
flow and that recombination takes place. For small clouds, 
the disruption takes a few 'cloud crushing' times, as de- 
fined by Klein, McKee and Colella (1994), by which time 
the leading shock will be a few times x^^"^ ahead of the 
cloud, where x is the density contrast between the cloud 
and intercloud medium. For a x < 100 this leads to a 
displacement of a few tens of the cloud size. As long as 
the density contrast or cloud sizes are not too large the 
optically emitting gas will be close to the boundary of the 
SNR and moving with a velocity a fraction of the initial 
shock speed. The difficulty with such a picture is that a 
very large number of small clouds are required to get the 
very smooth O ill distribution as seen in our figures. 

Note that the X-ray shock velocity estimate comes 
from the current X-ray temperature and if this has been 
cooled significantly by the addition of cold stationary 
clouds then the actual X-ray shock velocity should be even 
higher. However, McKee & Cowie (1975) argue that the 
acceleration and evaporation of the over-run clouds will 
have little effect on the energetics of the SNR as long as 
their filling factor is not too large. Some support for this 
model comes from the clumpy nature of the X-ray emis- 
sion as mapped by Lu & Aschenbach (2004). 

An alternative is that a pre-existing dense circumstel- 
lar shell has been encountered by the blast wave and this 
has led to the two different shock velocities. If the pro- 
genitor was a massive star then an H II region and wind 
blown shell will have been set up around the star. Shull 
et al. (1985) discussed the effects of a SN explosion in a 
cloiidy medium around a massive star. The formation of 
the H II region destroys any H I clouds and thus the SN 
explosion takes place in a cavity created by the star and 
then encounters a clumpy shell. Similarly, Charles, Kahn 
& McKee (1985) argue that the Cygnus loop SNR was 
caused by a SNR exploding in a pre-existing cavity, pre- 
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sumably generated by the progenitor. Such a model has 
the attraction that it accounts for the relatively low am- 
bient medium density implied by a Scdov analysis of the 
X-ray shock emission. It also offers a possible explanation 
for the morphology of G65.3+5.7, discussed further be- 
low, since an axisymmctric mass-loss rate prior to the SN 
explosion could later govern the evolution of the shape of 
the remnant. 

The local electron densities in the recombination zone, 
as measured in Paper 1, ranged between 30 — 170 cm~^. 
For a temperature of 10^ the thermal pressure is then 
nT ~ 3 X 10^ K cm-3 - 2 x 10^ K cm'^. In contrast, the 
thermal pressure of the X-ray gas implied by the shock 
speed and density measured by Lu fc Aschenbach (2004) 
is only 1.6 x 10^ K cm^'^. This is another X-ray/optical 
discrepancy that has been observed before: Raymond et 
al (1988) found a similar effect in the Cygnus loop SNR. 
There, the shock ram pressure exceeds the thermal pres- 
sure in both the optical and X-ray emitting gas. Raymond 
et al (1988) discuss a possibility that could be relevant 
here. If a blast wave impacts a large enough cloud, a re- 
verse shock is generated that halts the expansion of the 
rarefied X-ray emitting gas (McKee & Cowie 1975). This 
produces a temporary overpressure in the remnant that 
then accelerates the shock. The obstacle that causes this 
sequence of events could be a previously generated wind 
blown shell as discussed above. 

3.3. Morphology 

The more extended, deep image in Fig. lb reveals clearly 
the prominent, nearly circular, filamentary structure, « 
3° diameter, that must define the edge of the eastern lobe 
which is most probably expanding at w 155 km s^^. The 
prominent but less well-defined westerly lobe appears to 
be composed of at least two substructures and to com- 
plicate matters further a fainter northern lobe outside of 
the main perimeter of the remnant can now be seen. The 
northern perimeter of the easterly lobe is approaching the 
observer at around 16 km with respect to its southern 
counterpart. All that can be assumed is that the two edges 
have formed in separate interstellar clouds with different 
radial velocities. 

G65.3-I-5.7 appears to be a very good example of 
a barrel-shaped supernova remnant (see Kestevan & 
Caswell 1987 for a full discussion and many examples). 
There is a clear axis of symmetry, approximately east- 
west, either side of which are two bright limbs of emis- 
sion. Such a morphology is common and the mechanisms 
for producing it can be divided into extrinsic and intrin- 
sic effects (Gaensler 1998). An example of the former is 
a well ordered local ISM magnetic field that is then in- 
corporated into the shell of the SNR. An example of the 
latter is a previous axisymmetric outflow from the pro- 
genitor that then influences the evolution of the SNR (see 
above). Gaensler (1998) presents evidence that the bar- 
rel axes tend to be aligned with the galactic plane which 



then suggests an extrinsic effect (specifically, the galactic 
magnetic field stratifying the ambient medium) is respon- 
sible in many cases. The barrel axis of G65.3+5.7 is at 
approximately 45 degrees to the plane and the remnant 
is displaced from the plane by ~ 100 pc. This is of the 
order of the gas scale height and so the orientation of the 
field could be begining to depart from being parallel to 
the plane. Though the exact orientation is unclear, 45 de- 
grees is not unreasonable. It does seem unlikely however 
that interstellar clouds will be flattened and stratified at 
this angle to the ISM. One of the several proposed mech- 
anisms involving the compression of magnetic field (dis- 
cussed by Gaensler 1998) and the generation of the bright 
limbs could be responsible instead. 

4. Conclusions 

New images of G65.3-I-5.7 emphasize that this is a good 
example of a barrel-shaped supernova remnant. 

The high resolution line profiles presented here suggest 
that the global expansion is « 155 km s"-'^ which is around 
the shock velocity deduced from previous low-resolution 
spectra but well short of the 400 km s~^ forward shock 
velocity required to generate the X-ray emission. 

This expansion velocity combined with an expansion 
proper motion measurement (2.1 ± 0.4 arcsec in 48 yr) 
gives a distance of 770 ± 200 pc. 

It would be desirable to derive a more certain expan- 
sion velocity by obtaining further deep line profiles near 
the remnant's geometrical centre. The spectral resolution 
could be relaxed to « 30 km s^^ by broadening the slit 
width to ensure sufficient signal to noise ratio in a reason- 
able observing time for such faint nebulosity. 

Furthermore, a more accurate expansion proper mo- 
tion could be derived by obtaining current higher reso- 
lution (say 1") images of the filamentary edges though a 
filter that matches that of the baseline POSS image. 
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Fig. 1. A negative, continuum-subtracted mosaic of [O III] 5007 A images of G 65.3+5.7 is presented lightly in a) with 
the slit positions Pl-5 and 7-9 marked and where the prominent eastern and western lobes are apparent; and deeply 
in b) to reveal the fainter structure with a further northern lobe suggested (coords, are 2000 epoch). 
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Fig. 2. Negative greyscale presentations of the [O ill] 5007 A position-velocity arrays for EW slit positions (See Fig. 
la 

Pl-5 and 7-8. 
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Fig. 3. a) A small section of the image in Fig. 1 is shown with the full length of slit PI marked, b) The contour map 
is shown, with linear intervals, of a 43" section of the pv array for PI (see Fig. 2) where the slit intersects the bright 
filament. Here Bp iii] contours are separated separated by 1.25 x 10~^ erg s~^ cm~^ sr~^ A~^. 
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Fig. 5. [O III] 5007 A line profiles for the brightest filamentary regions intersected by the pv arrays in Fig. 2, and 
for the diffuse emission along slit P9, are presented. The profiles are for the various lengths (between the brackets) 
along the slits: for slit PI (7"), P2 (17"), P3 (61"), P5 (33"), P7 (104") and for P9 (113"- towards the westerly end 
of the slit position shown in Fig. la). The value given against each profile is for the peak surface brightness B[o iii] 
(uncorrected for interstellar extinction) in units of 10~® erg s~^ cm~-^ sr~^ A~^. 



